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ABSTRACT 

We derive the best characterization to date of the properties of radio quiescent massive galaxy 
clusters through a statistical analysis of their average synchrotron emissivity. We stacked 105 radio 
images of clusters from the 843 MHz SUMSS survey, all with Lx > 10^^ erg s~^ and redshifts z < 0.2, 
after removing point-source contamination and rescaling to a common physical size. Each stacked 
cluster individually shows no significant large-scale diffuse radio emission at current sensitivity levels. 
Stacking of sub-samples leads to the following results: (i) clusters with Lx > 3 x 10^^ erg s~^ show 
a 6cr detection of Mpc-scale diffuse emission with a 1.4 GHz luminosity of 2.4±0.4 xlO^^ W Hz-^ 
This is 1.5-2 times lower than the upper limits for radio quiescent clusters from the GMRT Radio 
Halo Survey (Venturi et al. 2008), and is the first independent confirmation of radio halo bi- modality; 
(n) clusters with low X-ray concentrations have a mean radio luminosity (2.6±0.6 xlO^^ W Hz~^) 
that is at least twice that of high X-ray concentration clusters, and (iii) both of these detections are 
likely close to the low- level "off-state" of GRHs in most or all luminous X-ray clusters, and not the 
contributions from a much smaller subset of "on-state" GRHs following the radio/X-ray luminosity 
correlation. Upcoming deep radio surveys will conclusively distinguish between these two options. We 
briefiy discuss possible origins for the "off-state" emission and its implications for magnetic fields in 
most or all luminous X-ray clusters. 
Subject headings: word - word - word 



1. INTRODUCTION 

Many massive clusters of galaxies host Mpc-scale dif- 
fuse radio emission that is not associated with ac- 
tive galactic nuclei (see IFerrari et al.l 120081 for a re- 
view). This synchrotron emission implies the pres- 
ence of juG magnetic fields, a population of cosmic-ray 
electrons/positrons (CRe^), and a spatially distributed 
mechanism to (re) accelerate the CRe^. The two most 
likely mechanisms for the maintenance of the CRe^ 
distribut ion are turbulent re-acceleration afte r a major 
merger (jBrunetti et al.l 120011 : iPetrosianI l2001l ) and sec- 
ondary production through collisions of longs:-live d cos- 
mic ray protons (CRp) with thermal protons (Den nisonI 
[l980"; 'Blasi & Colafrancesco'T999). The CRp can be ac- 
celerated at strong shocks during cluster formation, and 
resi de in the deep pot ential wells for a Hubble time (see, 
e.g., iBlasi et "al]l2007l for a review). 

Th e GMRT Radio Halo Survey (jVenturi et al.l I2007L 
l2008f ). which surveyed a complete sample of high X- 
ray luminosity clusters, showed a bi- modality in the 
GRH population which is consistent with expectations 
from models based on merger-driven re-acceleration 
(jBrunetti et al.ll2009f ). Those clusters that host a GRH 
follow a correlation between the radio and X-ray lumi- 
nosities, which we designate here as the "on-state" . The 
remaining "off-state" clusters show no signs of diffuse 
emission even in deep GMRT images. The distinguish- 
ing characteristic between on and off-state appears to be 
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the presence or absence of major m erger activity, respec- 
tively, although ou tliers do exist (|Cassano et al.l 120101 : 
iRussell et alll20TTI) . 

Despite the difficulties that pure secondary mod- 
els have in explaining recent observations (e.g., 
Brunetti etld] 2008: Donnert etal] l2010l : iBrown et all 
20111: iJelten^fe Profumo .201 there are reasons to 



expect Mpc-scale diffuse emission even in the off-state 
clusters. Though Fermi observations have placed upper 
limits on the ratio of non-thermal CRp to thermal en- 
ergy densities at t he level of ^5%-10% in several cl usters 
(Ackermann et all 120101 : iJeltema fc ProfumdIMTh . sec- 
ondary CRe^ from these allowed levels would still pro- 
duce faint diffuse emission in the presence of jaG mag- 
netic fields. Additionally, even off-state clusters in the 
turbulent re-acceleration model should host steep spec- 
trum GRHs that should be detectabl e with upcoming 
low frequency radio arravs fe.^.. iCas^no et al.l 120061 : 
iBrunetti et al.l 120081 : ICassano et al.l l2010bf). 

In this Letter, we provide the most stringent con- 
straints on the off-state synchrotron emission and GRH 
bi- modality by stacking images of massive, radio-quiet, 
clusters from the S ydney University Molonglo Sky Sur- 
vey (SUMSS, .Bock etld] [T999). SUMSS is a radio sur- 
vey of the Southern sky at 843 MHz, ~ 45'' resolu- 
tion, ~1 mJy beam~^ rms noise, and is very sensitive 
to extended low surface-brightness emission. We assume 
Ho = 70, = 0.7, = 0.3, and define the syn- 
chrotron spectral index a by 1^ ^ 

2. DATA ANALYSIS 

2.1. Initial Sample Selection 

We started with a sample of luminous X-ray clusters, 
Lx > 10^^ erg s~^, so that the clusters would be suffi- 
ciently massive to have deep potential wells, reservoirs of 
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thermal plasma, and have had time to accumulate CRp. 
We further restricted the clusters to z < 0.2; at that red- 
shift, the declination-dependent SUMSS resolution, with 
a maximum of 90'^ corresponded to ~295 kpc, allowing 
at least 3 independent beams across a 1 Mpc GRH. From 
the X -ray Galaxy Clusters Database (BAX, Sadat et aj] 
l2QQ4j) . there are 149 galaxy clusters with measured X- 
Ray luminosity Lx > 10^^ erg s~^, and with redshifts 
z < 0.2. 

2.2. Image Processing 

In order to uniformly probe diffuse Mpc-scale emis- 
sion through stacking, we first convolved each image as 
a function of redshift such that each beam is ~300 kpc 
in diameter. We then filtered out smaller-scale emission 
due to compact and slightly extended radio galaxies, us- 
ing the multi-scale spatial filter of IRudnickl (|2QQ2l) and 
a box size of 600 kpc. All images were re-scaled and 
resampled so that each pixel is the same physical size. 
After this processing, a 1 Mpc GRH of a given surface- 
brightness would appear the same at any redshift, and 
would add coherently if present in more than one (or 
all) of the clusters. Statistics for each image were calcu- 
lated (mean, rms, and peak fluxes) within a 1 Mpc box 
around the center, as well as a 1-3 Mpc annular box for 
background subtraction and noise estimation. 

2.3. Quality Control 

From these 149 processed images, we needed to elim- 
inate those that could have significant contamination 
from residual emission from bright point sources. To 
quantify the possible point source contamination, we in- 
jected synthetic 10 Jy sources into each of the images 
prior to the processing. The resulting stack showed 
~ 0.033 Jy beam~^ peak emission, or 0.33% leakage. 
Point source contamination in our stack is therefore neg- 
ligible compared to other noise contributions. 

Marginally extended sources could contaminate the fil- 
tered images at a higher level, so we eliminated their pos- 
sible influence empirically. After background removal, 
we plotted the diffuse, filtered flux within the 1 Mpc 
box against the peak unfiltered fiux within 1 Mpc (not 
shown). Clusters with peak fiuxes below 0.25 Jy beam~^ 
did not show any correlation between compact and fil- 
tered diffuse emission, and were kept for further analysis; 
105 clusters met this criterion. 

Our filtering procedure can, unfortunately, also reduce 
the observed fiux from our target diffuse emission. To 
assess the amount of missing fiux, we therefore injected 
1 Mpc, 0.5 Mpc, 0.25 Mpc, and 0.125 Mpc gaussians 
into each stacked image. We recover 91%, 60%, 18%, and 
4% of the total injected fiux, respectively. For the lumi- 
nosities listed in Table 1, we assumed 1 Mpc Gaussian 
halos, and thus scaled all the calculated luminosities up 
by 1/0.91=1.1. 

We also created a control sample of 105 randomly se- 
lected fields assigned to the same redshifts as the cluster 
sample, and processed those in an identical manner. 

Visual inspection of the cluster images did not reveal 
hints of large-scale diffuse emission. We compared in- 
tegrated fiux measurements from SUMSS and from the 
Parkes Telescop43 for several large sources, and found 

^ http:/ /www. parkes. atnf.csiro.au/ 
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Fig. 1. — Top: Measured radio luminosities for each cluster as a 
function of cluster X-ray luminosity. The majority are consistent 
with zero (dashed line), though a few show l-2cr positive and neg- 
ative deviations; Bottom: A plot of cluster Lx vs. c, with solid 
lines marking the regions of our sub-samples. 

that SUMSS is sensitive out to at least 20' scales (the 
closest cluster in our high-Lx sample is z=0.048, where 
1 Mpc is 17.9'). We further injected low signal-to- noise 
gaussian halos of various sizes directly into the images, 
and conclude that individual 1 Mpc halos with luminosi- 
ties of Pi. 4 GHz ~ 2 X 10^^ W/Hz would be visible in 
SUMSS over the range of redshifts in our sample (as- 
suming a=1.2). 

Figure [T] shows the distribution of the radio luminosi- 
ties (each individually a non-detection) as a function of 
Lx' Note that they can be either positive or negative due 
to the background subtraction. All radio luminosities 
(henceforth) are scaled to 1.4 GHz assuming an a = 1.2. 

2.4. Stacking & Subsamples 

Each image had somewhat different noise properties, 
as well as a redshift-dependent noise due to the spa- 
tial rescaling to a common physical size. We therefore 
stacked using a weighted mean, where the weight was the 
inverse variance (1/cr^) of each image in the background 
region between 1 and 3 Mpc. The results of stacking the 
105 cluster and control fields are shown in Figure [2] and 
Table 1. The cluster stack shows an ~4(j detection of 
diffuse emission of 1 Mpc in extent, with a peak fiux of 
0.99±0.1 mJy beam~^, while the control stack is consis- 
tent with noise. 

In order to compare with theoretical expectations, we 
then split the 105 clusters into subsamples based on two 
criteria; high vs. low Lx (> or < 3 x 10^^ erg s~^) , and 
level of X-ray concentration. 

To asses the concentration of X-rays we used the 
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parameter c ([Santos et al.l I2QQ8I : iCassano et al.l l2QlQf ), 
which is the ratio of the X-ray flux within r < 100 kpc to 
that within r < 500 kpc. Higher values of c indicate more 
centrally dominated, and thus relaxed, clusters, while 
lower values of c tend to indicate the presence of sub- 
structure induced by cluster mergers. We calcu lated c 
from the Rosat All Sky Survey (jVoges et alJll999f ) images 
of the clusters. We ordered the individual radio luminos- 
ity values of the 1043 clusters by their c parameter value 
and performed a Bayesian change-point analysis to deter- 
mine the value of c at which the distribution of luminosi- 
ties changed in mean and standard deviation. The anal- 
ysis found c = 0.156 as the most likely transition value, 
where the {mean, median, standard deviation) for clus- 
ters with c < 0.156 was (5.6, 2.4, 8.2) x 10^^ W Hz-\ and 
(2.1, 1.2, 26.0) X 10^^ W Hz-i for clusters with c > 0.156. 
We therefore use c = 0.156 to demarcate "merging" (low 
c) vs. "no n- merging" clusters. We note t hat this is very 
close to the value of c ~ 0.2 found by ICassano et al.l 
(|201Q) as the dividing line between halo and non-halo 
clusters. 

Figure [1] shows how the clusters are distributed in the 
Lx — c plane. The results from subsamples character- 
ized by either or both of Lx and c are shown in Figures 
[2] and Table 1. The high-Lx and low-c clusters both 
show a ~5<j detection of diffuse emission with 1 Mpc 
extent, while the \ow-Lx and high-c clusters are non- 
detections. Figure [3] shows the high-Lx luminosity plot- 
ted along with the GMRT RHS sample in the P - Lx 
plane, which is 1.5-2 times lower than the upper limits 
on radio quiescent clusters. 




TABLE 1 

Stack results for sub-sampling of the 105 cluster 

IMAGES. 



Lx [erg/s] 



PlA GHz [10"^ W/Hz] 



All 

> 3x10^4 

< SxlO"^"^ 

All 

All 



All 
All 
All 

<0.156 
>0.156 



105 
52 
53 
67 
37 



1.83±0.45 
2.37±0.39 
0.66±0.89 
2.62±0.66 
1.80±1.29 



^ Number of clusters in the stack 



3. DISCUSSION 

We first consider whether the origin of the diffuse emis- 
sion detected at high Lx is due to a small fraction of 
ordinary "on-correlation" radio halos that are just below 
the detection limit of the SUMSS survey, combined with 
a larger fraction of non-detections. To test this, we exam- 
ined the distribution of radio luminosity measurements 
(again, each individually a non-detection) for the high 
X-ray luminosity sample {Lx > 3 x 10^^ erg/s). We find 
that 10 of the 52 high-Lx clusters have radio errors con- 
siderably larger than the rest of the sample (see Fig. [1]), 
so we selected the 42 lowest noise clusters and their corre- 
sponding controls and compared their luminosity distri- 
butions. Application of the Kolmogorov-Smirnov (KS) 

^ one cluster was removed for this test due to high X-ray back- 
ground resulting in c > 1. 



Fig. 2. — For all images, greyscale runs from -1 mJy beam" 
(black) to +1 mJy beam"-*^ (white; beam= 300 kpc x 300 kpc), 
and contours start at 3a and increase at Icr intervals. See Table 1 
details of the samples. Images are 3 Mpc on a side. 

test gives a 5% probability that the clusters are drawn 
from the same population as the control sample (which 
by definition has zero contribution of real diffuse lumi- 
nosity), thus confirming the detection found in Fig. [2j 
The next question is whether the difference is due to a 
handful of luminous on-correlation clusters, or whether 
most / all clusters having a small positive signal far below 
the on-correlation. 

To address this, we first added a fixed luminosity to 
each control field and recalculated the probability that 
the cluster and control samples are drawn from the same 
population. The probability reaches a maximum of 91% 
when we add a luminosity of 5x10^^ W/Hz to each con- 
trol field, consistent with our measured stacked luminos- 
ity and the GMRT upper limits (Fig. [3j). 

We next performed a Monte Carlo analysis injecting an 
on-correlation luminosity into a small numbe r of the con- 
trol fie lds, using the radio halo fractions of iCassano et al.l 
(|2008l ). and leaving the rest unchanged. After 10^ ran- 
dom realisations, the mean likelihood that the observed 
and modified control distributions were the same was 
13±5%, with no cases >39%. While not statistically in- 
contestable, we do find it most likely that a large num- 
ber of clusters have an approximate luminosity of 3- 
5x10^^ W/Hz, representing their luminosity when they 
are in the "off-state". The fact that the average lumi- 
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Fig. 3. — Plot from lBrunetti et alJ ("2009") showing the correlation 
between GRH's radio and X-r ay luminosities. Arrows are upper 
limits from Venturi et al." (2008). The filled red circle is the stacked 
luminosity for the high-Lx sample, and the dashed red line is the 
value with the top three contributing clusters removed. Open red 
circles are results from stacking the top 13 and bottom 39 clusters 
in the high-Lx sample, all assuming 1 Mpc halos. The values are 
placed at the weight mean Lx-, where the weight is the predicted 
on-correlation value for each Lx (Brunetti et al. 2009). 

nosity of an off-state cluster is an order of magnitude 
less than on-correlation GRHs implies a bi-modal dis- 
tribution, consistent with results from the GMRT GRH 
Survey (e.g., Brunetti et al. 2009). 

We can further investigate how the stacked flux is dis- 
tributed among the clusters by plotting the cumulative 
luminosity. Figure 3] shows the cumulative luminosity of 
the entire cluster sample (black squares), the high-Lx 
sample (red circles), and the low-c sample (blue crosses), 
all of which produced positive detections (Fig. [2]). The 
cumulative luminosity was calculated in ascending Lx 
order, and we have applied the same weighting as in our 
stacking. There are several important features of the 
distribution: 1) There are several locations where the 
distribution jumps, indicating strong contributions from 
a few clusters; 2) In addition to these jumps, the over- 
all positive trend above Lx > 3 x 10^^ erg/s indicates 
a significant distribution of the flux among most of the 
clusters (as demonstrated statistically above); 3) The off- 
state emission is a function of Lx, even within the crude 
sub-sampling performed in the stacks. 

The Lx dependence can be investigated by further 
splitting the high-Lx sample into two. The top 13 clus- 
ters show a strong linear rise in cumulative luminosity, 
so we stacked the top 13 and the remaining 39 sepa- 
rately, and show the images in Fig. |4] and values in Fig. 
[3l We see significant emission from both subsamples, 
again confirming that the off-state emission comes from 
many clusters. We note that the value obtained in the 
39 cluster stack suggests that bi-modality extends down 
to lower Lx clusters than has been probed by previous 
studies (see Fig. [3]). 

The statistics of GRHs observed with the GMRT 
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Fig. 4. — Plot of the cumulative luminosity, noise- weighted as in 
our stacks, of the entire sample (black boxes), high-Lx clusters (red 
circles), and low-c clusters (blue crosses). Shaded regions indicate 
the sub-stacking of the top 13, middle 39, and bottom 53 clusters 
in the sample (see text), and the stacked images are shown as insets 
with the same greyscale and contours as Fig. [2] 

suggest that, if our sample where complete, we would 
expect roughly 2-4 sources in the high-Lx sample 
(jCassano et al.|[2QQ8f ). We therefore took a more conser- 
vative approach for calculating the off-state and assumed 
that there could be as many as 3 GRHs in our sample. 
We removed the top three contributors to the cumula- 
tive distribution and recalculated the stack, obtaining 
PiAGHz = 1.95 ± 0.75 X 10^^ W/Hz. We believe this 
value is a better estimate of the true off-state luminosity 
for high-Lx clusters. 

The presence of low-level Mpc-scale emission in 
off-state clusters requires a cluster-wide CRe^ acceler- 
ation mechanism and minimum ener gy ICM magnetic 
fields of ~ 1 /iG (assuming a = 1.2; iGovoni fc Ferettil 
l2QQ4l ). A similar claim of ubiquitous cluster magnetic 
fields is made in several studies based on ro tation- 
measure/depolarization (e.g., Ca rilli fc Taylorl I2QQ2I : 
Clarke 2004; Bonafede et al. 2011), although the ad- 
equate size samples of background sources needed to 
conduct a fair test are still missing. Possible origins for 
the CRe^ include cluster-wide acceleration processes 
such as turbulent re-acceleration and CRp secondaries, 
or emission from diffuse radio galaxies. We consider 
each of these below. 

Turbulent Re- acceleration: There is a large body 
of evidence connecting the presence of GRHs wit h cluster 
mergers (Buote 2001; Govoni et al. 2004; Cassan Q"et al.l 
120101) . Our results (Fig[2j) show that disturbed clusters 
have stronger off-state Mpc-scale emission. However, 
Fig. [T] indicates that this is partially de g enera te with 
high-Lx for our sample. iCassano" etall (|20T0h found 
that on-state radio halos were not found in clusters 
with c > 0.2, the radio luminosity of these systems 
being constrained by upper limits from GMRT radio 
observations at 610 MHz.. Prima facie^ our results show 
that high Lx and disturbed clusters preferentially show 
Mpc-scale diffuse emission in off-state clusters as well. 
It has been suggested that radio halos with very steep 
spectrum are generated in connection with less energetic 
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merger events (e.g., iCassano et al.|[2QQ6l : iBrunetti et al.l 
I2OO8). These sources should appear preferentially at 
lower radio frequencies and would be sub-luminous halos 
at ~GHz frequencies. Thus one possibility is that our 
stacked radio signal from off-state clusters is from these 
faint systems. LO FAR is expected to detect most of 
these systems (e.g.. ICassano et al.ll2QlQlJ ) and will allow 
for testing of this possibility. 

CR-PROTONS: An exciting possibility is that we are de- 
tecting diffuse emission from secondary CRe^ resulting 
from CRp collisions w ith thermal protons. Fermi limits 
(jAckermann et an i2QlQ) allow the CRp energy density 
to be at most 5-10% of the thermal energy density in 
clusters, while GMRT upper limits from the-non detec- 
tion of diffuse Mpc-scale emission in off-state clusters 
constrain the ratio of CRp and thermal energy density 
to less than several percent, provi ded these systems ar e 
magnetized at the few /iC level (jBrunetti et al.ll2QQ7[ ). 
The level of emission detected in these systems is close to 
these GMRT limits, implying a ratio of CRp to thermal 
energy density ~3% (in 1 Mpc^ volume), assuming 
< B >^ 2 /iC. In a recent paper Brunetti & Lazarian 
(|2Qlll ) calculated models where primary CRp and their 
secondary products are reaccelerated by merger-driven 
MHD turbulence generating giant radio halos. Accord- 
ing to these calculations the synchrotron emission from 
secondary particles in off-state non-turbulent systems 
is a factor ~10 fainter than that in GRH, in rough 
agreement with our observations. In this scenario it 
is also expected that the synchrotron emission scales 
with the amount of turbulence in the hosting clusters, 
with GRH becoming gradually more powerful (and with 
harder spectra) in more turbulent (disturbed) systems. 
This may explain why we detect Mpc-scale emission 
preferentially in low-c clusters. Additional possibilities 
for understanding the correlation with cluster dynamics 
might stem from differences in magnetic properties 
and CRp dy namics in relaxed and disturbed clusters 
(|Keshet fc Lo eb 2010; Keshet 2010; EnBlin et all Uoil). 

Diffuse Radio Galaxies: It is possible that these 
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massive clusters host faint synchrotron emission from 
past/current AGN activity, though we don't see any ac- 
tively driven jets on these scales in the individual images. 
Absent the jet driven power, the CRe^ would still re- 
quire some in situ re-acceleration in order to be present 
on these large scales. Additionally, the emission is corre- 
lated with disturbed X-ray morphology (low-c), while it 
is typically cooling-core clusters (high-c) that host pow- 
erful extended radio galaxies in their core. These sources 
are also heavily suppressed by our filtering. 

4. SUMMARY 

We performed a statistical analysis of the diffuse, Mpc- 
scale, radio emission in 105 massive clusters from the 
SUMSS survey. Our major findings are 1) Mpc-scale 
diffuse radio emission is detected in formerly radio qui- 
escent clusters, at a level slightly below current upper 
limits. This emission is likely originating from clusters in 
the off-state; 2) there is a substantial population of clus- 
ters with Pi. 4 GHz much less than the observed P — Lx 
luminosity correlation for GRHs, confirming GRH bi- 
modality. We further confirm bi-modality in lower Lx 
clusters than has previously been known; 3) off-state 
emission is correlated with high-Lx and di sturbed cluster 
morphology, analogous to on-state GRHs (jCassano et al.l 
l2QlQf ): 4) Magnetic fields at the ~1 /iG level are present 
on Mpc-scales even in formerly radio-quiescent massive 
clusters. We demonstrate the power of cluster stacking, 
which will allow us to probe lower mass/Lx clusters and 
groups wit h upcoming all-sky radio continuum surveys 
(e.g., Norris et al.ll2Qll[ ). 
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